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We have searched for the decays —> Dfa^, B^ —> Dfaif and B^ in 

a sample of about 230 million T(4S)—>-B-B decays collected with the BABAR detector at the PEP-II 
asymmetric-energy B Factory at SLAC. We find no evidence for these decays and set upper limits at 
90% C.L. on the branching fractions: B{B° Dfa^) < 1.9 x 10“®, B{B^ —> D^'^a^) < 3.6 x 10~®, 
B{B° D+a^) < 1.9 X lO"**, and B{B° D*+a^) < 2.0 x 10“'‘. 


PACS numbers: 13.25.Hw, 12.15.Hh, ll.30.Er 


The time-dependent decay rates for neutral B mesons 
into a D meson and a light meson provide sensitivity 
to the Cabibbo-Kobayashi-Maskawa (CKM) Q quark 
mixing matrix phases (3 and 7 |3- A CP-violating 
term emerges through the interference between 
mixing mediated and direct decay amplitudes. The 
time-dependent CP-asymmetries in the decay modes 
P° ^ 13 have been studied by BABAR and 

BELLE Hi In these modes, the CP-asymmetries arise 
due to a phase difference between two amplitudes of very 
different magnitudes: one decay amplitude is suppressed 
by the product of two small CKM elements Vub and Vcd, 
while the other is CKM favored. Therefore, the decay 
rate is dominated by the CKM-favored part of the ampli¬ 
tude, resulting in a very small CP-violating asymmetry. 

Recently it was proposed to consider other types of 
light mesons in the two-body final states [3 • The idea is 
that decay amplitudes with light scalar or tensor mesons, 
such as a(}" or aj, emitted from a weak current, are sig- 
nihcantly suppressed because of the small coupling con¬ 
stants fao( 2 )- the SU{2) limit, fa^ = 0 (since the 
coupling constant of a light scalar is proportional to the 
mass difference between u and d quarks), and any non¬ 
zero value of fag is of the order of isospin conservation 
breaking effects. Since the light tensor meson has spin 
2, it cannot be emitted by a VE-boson (i.e. fa 2 = 0), and 
thus could only appear in a ICfc-mediated process via final 
state hadronic interactions and rescattering. Therefore, 
the absolute values of the CKM-suppressed and favored 
parts of the decay amplitude (see Figure Q top two di¬ 
agrams) could become comparable, potentially resulting 
in a large CP-asymmetry. No B —>■ ao( 2 )-^ transitions 
have been observed yet. A summary of the theoretical 
predictions for the values of Vub and Kcfc-niediated parts 
of the P° ^ branching fractions can be found 

in 0. 

The Tufa-mediated amplitudes in [3 were computed in 
the factorization framework. In addition to model uncer¬ 
tainties, signihcant uncertainty in the theoretical calcu¬ 
lations is due to unknown B ao( 2 )^ transition form 
factors. One way to verify the numerical assumptions 
and test the validity of the factorization approach ex¬ 


perimentally is to measure the branching fractions for 
the SU{3) conjugated decay modes P° ^ Pi*^^ao( 2 )- 
These decays are represented by a single tree diagram 
(Figure n bottom diagram) with external VE+ emission, 
without contributions from additional tree or penguin di¬ 
agrams. The t4fa-mediated part of the P° ^ 

decay amplitude can be related to P° ^ us¬ 

ing tan (6>cabibbo) = \Vcd/Vcs\ and the ratio of the decay 
constants fjjM/foM- 

Branching fractions of P° —*■ are predicted 

to be in the range 1.3-1.8 (2.1-2.9) in units of 10“® 13 - 
Branching fraction estimates for P° ^ of ap¬ 

proximately 8 x 10“^ are obtained using SU (3) symmetry 
from the predictions made for P° ^ D^*'>~^aQ in |3|. 



FIG. 1: Top diagrams: tree diagrams contributing to the 
decay amplitude of B^ —> cIq(^ 2 ) (including the B°B^ 

mixing mediated part of the amplitude). Bottom diagram: 
tree diagram representing the decay amplitude of B° 

1^3 “o( 2 )- 

In this paper we present the hrst search for the de¬ 
cays P° ^ P+Oq, ^ P*+aQ, ^ and 

B^ —>• The analysis uses a sample of approxi¬ 

mately 210 fb“^, which corresponds to about 230 mil¬ 
lion 7’(4S') decays into BB pairs collected in the years 
1999-2004 with the BABAR detector at the asymmetric- 
energy P-factory PEP-II |3] ■ The BABAR detector is de¬ 
scribed elsewhere and only the components crucial 
to this analysis are summarized here. Charged particle 
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tracking is provided by a five-layer silicon vertex tracker 
(SVT) and a 40-layer drift chamber (DCH). For charged- 
particle identification, ionization energy loss (dE/dx) in 
the DCH and SVT, and Cherenkov radiation detected 
in a ring-imaging device are used. Photons are identi¬ 
fied and measured using the electromagnetic calorimeter, 
which is comprised of 6580 thallium-doped Csl crystals. 
These systems are located inside a 1.5 T solenoidal su¬ 
perconducting magnet. We use GEANT4 EH software 
to simulate interactions of particles traversing the BABAR 
detector, taking into account the varying detector condi¬ 
tions and beam backgrounds. 

The selection criteria are optimized by maximizing the 
ratio of expected signal events S to the square-root of the 
sum of signal and background events B. For the calcula¬ 
tion of S we assume B{B^ to be the mean 

values of the predicted intervals from [Sj and an estimate 
of B{B'^ ^ obtained from B{B^ ) 

predicted in [2| and assuming SU (3) symmetry. The op¬ 
timal selection criteria as well as the shapes of the distri¬ 
butions of selection variables are determined from sim¬ 
ulated Monte Carlo (MC) events. We use MC samples 
of our signal modes and, to simulate background, inclu¬ 
sive samples of B^B~ (800 fb“^), B°B^ (782 fb“^), cc 
(263 fb“^), and gg, q = u,d,s (279 fb“^). In addition, we 
use large samples of simulated events of rare background 
modes which have final states similar to the signal. 

Candidates for Df mesons are reconstructed in the 
modes Df —> (/)7r+, and with (j) 

K~^K~, K*^ K~'k~^ and > tt+tt”. The can¬ 

didates are reconstructed from two oppositely-charged 
tracks , w ith an invariant mass close to the nominal 
mass |l2j. that come from a common vertex displaced 
from the e“'"e“ interaction point. All other tracks are re¬ 
quired to originate less than 1.5 cm away from the e“'"e“ 
interaction point in the transverse plane and less than 
10 cm along the beam axis. Charged kaon candidates 
must satisfy kaon identification criteria that are typically 
around 95% efficient, depending on momentum and po¬ 
lar angle, and have a misidentification rate at the 10% 
level. The cj) K'^K~, K*^ K~Tr^ and tt+tt" 

candidates are required to have invariant masses close to 
their nominal masses (we require the absolute dif¬ 
ferences between their measured masses and the nominal 
values to be in the range 12-15 MeV, 35-60 MeV 
and 7-12 MeV, respectively, depending on the and 
Df decay modes). The polarizations of the K*^ and (j) 
mesons in the Df decays are used to reject backgrounds 
through the use of the helicity angle 9h, defined as the 
angle between the K~ momentum vector and the direc¬ 
tion of flight of the Df in the K*° or cj) rest frame. The 
K*° candidates are required to have | cos 9h \ greater than 
0.25-0.5 and (f) candidates are required to have | cos9h\ 
greater than 0.3-0.5, depending on the decay mode. 
We also apply a vertex fit to the Df candidates that 


decay into (^7r+ and A'*°Ar+, since all charged daughter 
tracks of Df are supposed to come from a common ver¬ 
tex. The of the vertex fit is required to be less than 
10-16 (which corresponds to a probability of better than 
0.1% — 1.9% for the 3 track vertex fit), depending on the 
reconstructed mode. 

The candidates are reconstructed in the mode 
> D+y. The photons are required to have an 
energy greater than 100 MeV. The Df and D*~^ can¬ 
didates are required to have invariant masses less than 
about ±2(7 from their nominal values 0. The invariant 
mass of the is calculated after the mass constraint 
on the daughter Df has been applied. Subsequently, all 
candidates are subjected to a mass-constrained fit. 

We reconstruct and candidates in their decay 
to the r]TT~ final state. For reconstructed ry ^ 77 candi¬ 
dates we require the energy of each photon to be greater 
than 250 MeV for Oq candidates, and greater than 300 - 
400 MeV for candidates, depending on the Df mode. 
The r] mass is required to be within a ± 1(7 or ± 2 t 7 interval 
of the nominal value na , depending on the background 
conditions in a particular i?°, Df decay mode (the 77 
mass resolution is measured to be around 15 MeV/c^). 
The a)}' and aj candidates are required to have a mass 
in the range 0.9-1.1 GeV/c^ and 1.2-1.5 GeV/c^, 
respectively. We also require that photons from 77 and 
D*'^ are inconsistent with hypothesis when combined 
with any other photon in the event (the tt^ veto window 
varies from ±10 to ±15 MeV/c^). Finally, the B^ meson 
candidates are formed using the reconstructed combina¬ 
tions of D'^Qq, Dg+ag and 

The background from continuum gg production (where 
q = u,d, s, c) is suppressed based on the event topology. 
We calculate the angle {9t) between the thrust axis of 
the B meson candidate and the thrust axis of all other 
particles in the event. In the center-of-mass frame (c.m.), 
BB pairs are produced approximately at rest and have 
a uniform cos 9 t distribution. In contrast, gg pairs are 
produced in the c.m. frame with high momentum, which 
results in a | cos0t| distribution peaking at 1. Depending 
on the background level of each mode, | cos6>t| is required 
to be smaller than 0.70-0.75. We further suppress back¬ 
grounds using a Fisher discriminant (E) 0 constructed 
from the scalar sum of the c.m. momenta of all tracks 
and photons (excluding the B candidate decay products) 
flowing into 9 concentric cones centered on the thrust 
axis of the B candidate. The more isotropic the event, 
the larger the value of T. We require T to be larger than 
a threshold that retains 75% to 86% of the signal while 
rejecting 78% to 65% of the background, depending on 
the background level. In addition, the ratio of the second 
and zeroth order Fox-Wolfram moments 0 must be less 
than a threshold in the range 0.25-0.40 depending on the 
decay mode. 

We extract the signal using the kinematical variables 
- (E* and AA = Xi* \/wf ± pf - 
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where is the beam energy in the c.m. frame, p* is 
the c.m. momentum of the daughter particle i of the 
meson candidate, and is the mass hypothesis for parti¬ 
cle i. For signal events, m-ES peaks at the meson mass 
with a resolution of about 2.7 MeV/c^ and AE peaks 
near zero with a resolution of 20 MeV, indicating that 
the B^ candidate has a total energy consistent with the 
beam energy in the c.m. frame. The B^ candidates are 
required to have | Ai?| < 40 MeV and wes > 5.2 GeV/c^. 

The fraction of multiple B^ candidates per event is es¬ 
timated using the MC simulation and found to be around 
2% for Tl^aQ( 2 ) and 5% for combinations. In 

each event with more than one B^ candidate that passed 
the selection requirements, we select the one with the 
lowest |Ai?| value. 

After all selection criteria are applied, we estimate the 
B^ reconstruction efficiencies, excluding the intermediate 
branching fractions (see Table P). 


TABLE I: Reconstruction efficiencies for 5° —> 
decays (exclnding the intermediate branching fractions). 


Decay mode 

P/ ^ 

P+ ^ K*°K+ 

P+ ^ K°K+ 

B° ^ P+ao 

4.7% 

2.9% 

2.5% 

B° D+a^ 

1.9% 

1.1% 

1.1% 

^ Dt+a^ 

2.2% 

1.5% 

1.3% 

^ D*+a2 

0.9% 

0.7% 

0.5% 


Background events that pass these selection criteria are 
mostly from qq continuum, and their tues distribution is 
described by a threshold function [l^ : 

/(toes) ~ Toes vl-^exp [-5(1 - 

where x = 2toes / y/s, y/s is the total energy of the beams 
in their center of mass frame, and 5 is the fit parame¬ 
ter. A study using simulated events of B^ and decay 
modes with final states similar to our signal mode, in- 
eluding Ds tt~ and Dl ’ p“, shows that these modes 
do not peak in toes- 

Figure [3 shows the toes distributions for the recon¬ 
structed candidates 5° ^ D+Og , B^ 5° ^ 

D*'^Qq and ^ D*'^a 2 - For each mode, we perform 
an unbinned maximum-likelihood fit to the toes distri¬ 
butions using the candidates from all Df decay modes 
combined. We ht the toes distributions with the sum 
of the function f{mEs) characterizing the combinatorial 
background and a Gaussian function to describe the sig¬ 
nal. The total signal yield in each B^ decay mode is 
calculated as a sum over modes {i = (('7r+, 

K°K+y. 

^sig ~ ^ ^BB ' ^ ^ ‘ 


where B is the branching fraction of the B^ decay mode, 
EEE is the number of produced BB pairs, Bi is the prod¬ 
uct of the intermediate branching ratios and is the re¬ 
construction efficiency. The mean and the width of the 
Gaussian function are fixed to values obtained from sim¬ 
ulated signal events for each decay mode. The threshold 
shape parameter 5, along with the branching ratio B are 
free parameters of the fit. The likelihood function is given 
by: 


-N N 

^ + iN- ®), 


where P/*® and are the probability density functions 
for the corresponding hypotheses, N is the total number 
of events in the ht and i is the index over all events in 
the ht. 



mES (GeV/c^) 

FIG. 2: Distributions of toes for B° —> Di*P candi¬ 
dates overlaid with the projection of the maximum likelihood 
fit. Contributions from P/ modes are shown with a different 
hatching style. The fit procedure and results are described in 
the text. 


Table m (second column) shows the signal event yields 
from the toes fit. Due to a lack of entries in the signal 
region for the B° —> P*+a 2 mode, the ht did not yield 
any central value for the number of signal events in this 
mode. Accounting for the estimated reconstruction ef- 
hciencies and daughter particles branching fractions, we 
measure the branching fractions shown in the third col¬ 
umn of Table QH 

The systematic errors include a 14% relative uncer¬ 
tainty for Df decay rates ^ . Uncertainties in the toes 
signal and background shapes result in 11% relative er¬ 
ror in the measured branching fractions. The rest of the 
systematic error sources, which include uncertainties in 
photon and ry reconstruction efficiencies, the a^}" and 
masses and widths, track and AT® reconstruction, charged 
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TABLE II: Signal yields, branching fractions and upper lim¬ 
its on the branching fractions for decays. 

Numbers in parentheses in the third and fourth columns indi¬ 
cate the branching fractions and the upper limits multiplied 
by the branching fractions of the decays Df —> and 

4(2) ^ 


B° mode 

'^sig 

B [10-®(10-'^)] 

U.L. [10-®] 

U+tto 


0.6t);t±0.1 (2.61®;® ±0.5) 

1.9 (0.09) 

D+a^ 

U.O_0.6 

6Atl°r^±l.5 (4.5l4;g ±0.8) 

19 (0.13) 

D ^ <1q 

1 c:+2.3 
-*-•^-1.8 

1.4t2-i±o.3 (6.5t7°g^ ± 1.2) 

3.6 (0.17) 

D*+a^ 

- 

- (-) 

20 (0.13) 


kaon identification, range between 3% and 10%. We as¬ 
sume the branching fraction for aj ^ 7771 + to be 100 % 
and assign an asymmetric systematic error of — 10 % to 
this assumption. The systematic error in the number 
of produced BB pairs is 1.1%. It was checked that the 
selection of the best candidate based on |Ai?| does not 
introduce any signihcant bias in the toes fit- The to¬ 
tal relative systematic errors are estimated to be around 
25% for each mode. 

We use a Bayesian approach with a flat prior above 
zero to set 90% confidence level upper limits on the 
branching fractions. In a given mode, the upper limit 
on the branching fraction (Bul) is defined by: 

rBuL P<x> 

/ C{B)dB = 0.9 X / £{B)dB 

Jo Jo 

where C{B) is the likelihood as a function of the branch¬ 
ing fraction B as determined from the toes fit described 
above. We account for systematic uncertainties by nu¬ 
merically convolving C{B) with a Gaussian distribution 
with a width determined by the relative systematic un¬ 
certainty multiplied by the branching fraction obtained 
from the toes fit. In cases with asymmetric errors we 
took the larger for the width of this Gaussian function. 
In case of (where no central value was determined 

from the ht) we conservatively estimate the absolute sys¬ 
tematic error by taking the numerically calculated 90% 
confidence level upper limit (without the systematic un¬ 
certainties) instead of the fitted branching fraction. The 
resulting upper limits are summarized in Table ITTI (fourth 
column). The likelihood curves are shown in Figure 0 

We have also calculated upper limits without including 
the intermediate branching fractions of the decays D'^ 

0 and a)|) 2 ) ^ E3- The relative systematic 

errors in this case are reduced to 18% for each of the 5° 
meson decay modes. The results are presented in Tablellll 
(third and fourth columns, numbers in parenthesis). 

In conclusion, we do not observe any evidence for 
the decays B^ ZD+Oq, 5° ^ Df ^ II*+ag 



FIG. 3: Likelihood functions of the fit for the toes distri¬ 
butions of the selected candidates. Solid 

curves represent the original likelihood scan from the fit, 
the dashed lines show the result of the convolution with 
the systematic errors Gaussian. Vertical lines indicate the 
90% Bayesian C.L. upper limit value. 


and B^ —> D*~^a 2 , and set 90% G.L. upper limits on 
their branching fractions. The upper limit value for 
5° ^ Dt% is lower than the theoretical expectation, 
which might indicate the need to revisit the B —> gqX 
transition form factor estimate. It might also imply the 
limited applicability of the factorization approach for this 
decay mode. The upper limits suggest that the branch¬ 
ing ratios of B^ are too small for CP- 

asymmetry measurements given the present statistics of 
the i?-factories. 
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